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<^ ■ Abstract 

^^ . Contractile forces are essential for many developmental processes involving cell shape change and 

tissue deformation. Recent experiments on reconstituted actomyosin networks, the major compo- 
nent of the contractile machinery, have shown that active contractility occurs above a threshold 
O ■ motor concentration and within a window of crosslink concentration. We present a microscopic 

^N^, dynamic model that incorporates two essential aspects of actomyosin self-organization: the asym- 

Qh" metric load response of individual actin filaments and the correlated motor-driven events mimicking 

myosin-induced filament sliding. Using computer simulations we examine how the concentration 
\^ , and susceptibility of motors contribute to their collective behavior and interplay with the network 

\Q ' connectivity to regulate macroscopic contractility. Our model is shown to capture the forma- 

^ ■ tion and dynamics of contractile structures and agree with the observed dependence of active 

^vq . contractility on microscopic parameters including the contractility onset. Cooperative action of 

load-resisting motors in a force-percolating structure integrates local contraction/buckling events 
K> ■ into a global contractile state via an active coarsening process, in contrast to the flow transition 

C^ ■ driven by uncorrelated kicks of susceptible motors. 



I. INTRODUCTION 

Contractile forces are essential for many processes vital to development, ranging from 
cytokinesis and cell motility^ to wound healing and gastrulation^. Networks of filamentous 
actin (F-actin) and the molecular motor, type II myosin have been identified as the major 
components of the contractile machinery. The actin network provides a structural scaffold on 
which the myosin motors move, powered by ATP hydrolysis. Actomyosin networks generate 
contractile forces through the activity of myosin motors, which themselves assemble into 
bipolar minifilaments that generate sustained sliding of neighboring actin filaments relative 
to each other in order to reorganize F-actin networks and generate tension.- When coupled 
to the cell substrate or using cell-cell adhesions, contractile actomyosin networks transmit 
forces to their environment. 

In addition to the microtubule-kinesin system, another important filament-motor as- 
sembly in cells that forms well-focused mitotic spindle poles driven by polarity sorting 
mechanism^^^ to accomplish high-accuracy segregation of duplicated chromosomes, acto- 
myosin condensates appear in diverse tissues and organisms as transient structures that 
coalesce into still larger arrays that exert contractile forces. Examples include the contrac- 
tile rings driving cytokinesis and wound healing, and the contractile networks that deform 
epithelial cell layers in developing embryos and drive polarizing cortical fiows^*^. 

Some recent theoretical efforts have modeled the contractile actin cortex as an active po- 
lar gel and have derived effective continuum theories within a hydrodynamic framework.-"— 
These macroscopic approaches based on generic symmetry considerations predict the for- 
mation of diverse patterns in acto- myosin gels such as asters and ring-like structures, which 
have been observed in studies in vitro.— 

Recently Bendix and coworkers have reconstituted contractility in a simplified system 
of F-actin, muscle myosin II motors, and a-actinin crosslinks.— The well- controlled na- 
ture of this in vitro system allows a systematic study of the dependence of contractility 
on microscopic parameters, such as the number and activity of myosin motors, crosslink 
density and actin network connectivity. It has been shown that contractility occurs above 
a threshold motor concentration and within a window of crosslink concentrations. Whereas 
earlier experiments on purified actomyosin solutions have established that contraction of 
F-actin networks by myosin II motors at physiological ATP level requires the presence of F- 



actin crosslinks,— li^ as has been confirmed by recent theoretical work,^^ this newly observed 
non-monotonic dependence of contraction tendency on crosslinking strength still calls for 
explanation. 

There are two important aspects in actomyosin self-organization: (1) Actin filaments have 
a highly asymmetric response to axial loading: They strongly resist tensile forces but easily 
buckle under compressive loads of several piconewtons. The ability to sustain large tension 
allows the motor-induced stresses to propagate significant distances through the network, 
whereas the buckling instability promotes formation of local actomyosin aggregates that 
coalesce into larger arrays exerting contractile forces (as exemplified in the spontaneous for- 
mation of myosin foci starting from a uniform distribution in an isotropic actin network—*^). 
(2) Motor-induced movements come in correlated pairs: the motor acts on a pair of parallel 
filaments to slide them past one another, inducing a pair of equal and oppositely directed 
moves at the crosslinks which are thus pulled together. Despite the success of hydrody- 
namic theories in predicting diverse patterns, a more microscopic model that could capture 
the nonlinear buckling behavior and the correlated motor-driven events, both of which are 
crucial for actomyosin self-organization and active contractility, is needed. 

Here we provide a microscopic dynamic model for active contractility that combines the 
motor-driven stochastic processes (modeled as correlated kicks on motor-bonded crosslinks) 
with the asymmetric load response of individual actin filaments. This minimal model will 
be shown to exhibit the experimentally observed dependence of macroscopic contraction on 
motor concentration and actin network connectivity. The model also highlights the key role, 
in structural development, of motor susceptibility, a parameter characterizing how sensitively 
the motors respond to imposed forces. 

By performing dynamic Monte Carlo simulations, we investigate the formation and 
dynamics of nonequilibrium structures in an actomyosin network modeled as a "cat's 
cradle'—*^ consisting of crosslinked nonlinear-elastic filaments subject to anti-correlated 
kicks on motor-bonded crosslinks (see Fig. [1] for a schematic illustration). We first study 
how the concentration and susceptibility of motors determine the collective behavior gen- 
erating diverse patterns. We further construct a phase diagram for active contractility as 
a function of motor concentration and network connectivity at a given motor susceptibil- 
ity. This diagram identifies the threshold motor concentrations and contains a window of 
network connectivity for active macroscopic contraction, consistent with observation.— We 



also find that at high connectivity contraction can still occur (at intermediate motor con- 
centrations) but only if the excluded volume effect is negligibly small. We finally compare 
the structures that develop for systems with correlated kicks to those when there are only 
uncorrelated kicks as in our earlier work^^. In particular, under uphill-prone motor kicks, 
the formation of "asters" is replaced by the formation of disordered condensates resulting 
from an active multistage aggregation process due to buckling of connected actin structures 
induced by cooperative anti-correlated kicks. This prediction receives support from recent 
in vitro experimental studies on how the collective action of myosin motors organizes actin 
filaments into contractile structures.— 

Multiple factors that cause dynamic remodeling in cells are clearly absent in the re- 
constituted assay,— such as the disassembly of contractile structures and the transience of 
physiological actin cross-linking proteins. In our numerical model, we likewise assume that 
network connectivity and motor distribution are quenched once initially assigned, in line 
with the fact that in vitro structures are irreversibly assembled because many in vivo fac- 
tors allowing fast pattern renewals are left out. Therefore the nonequilibrium dynamics and 
structures exhibited in our numerical study arise solely from the intrinsic activity of motors 
firmly built into the actin network driving correlated movements stochastically. This con- 
trasts with the two-fluid model that treats the cytoskeletal network as an elastic continuum 
where motor (un) binding kinetics leads to enhanced low- frequency stress fluctuations.— The 
discreteness of the power strokes and thus the kick steps in our model bears relevance to the 
pulsed contraction observed in actomyosin networks in vivo.— 



II. MODEL 



A. Model system and dynamic rule 



We model the actomyosin network as a cat's cradle^ii^^ consisting of nonlinear elastic 
filaments built on a three dimensional random lattice of volumeless crosslinks (Fig. [1]). 
These nonlinear filaments stretch elastically with effective stiffness /37 when their con- 
tour length r exceeds the relaxed length Lg but buckle and become fioppy upon short- 
ening, as described by the pair interaction potential between bonded neighboring crosslinks 
(5U{r) = B(r — Le)l3^{r — I/e)^/2 where /3 = l/fc^T and B(-) is the Heaviside step func- 
tion. The assumed weakness of the excluded volume effect allows large-scale structural 
rearrangements. We assign two mean field parameters to characterize the architecture of 
this filament-motor assembly: (1) the network connectivity Pc, which denotes the fraction 
of nearest-neighbor pairs of crosslinks bonded by filaments; (2) the motor concentration, P^? 
which indicates the fraction of active bonds, i.e. those attached by motors which induce 
equal and oppositely directed kicks on the connected crosslink pair. Network connectivity 
and motor concentration determine the fraction of active nodes against passive nodes: active 
nodes have motor- attached bonds and are subject to anti-correlated kicks with their motor- 
bonded neighbors whereas passive nodes have no motor-attached bonds and only undergo 
Brownian motion. 

To mimic the motor-driven filament sliding in actomyosin networks we describe the mo- 
tors as generating anti-correlated kicks on pairs of crosslinks along their lines of centers. 
Assuming a fixed kick step size /, consistent with the nearly periodic structure of the 
actin filaments, an anti-correlated kick pair acting on nodes i and j can be represented 
by (/^j, Iji) = l{fij^ ~fij)) where fij is a unit vector pointing from node i to node j. These 
anti-correlated kick pairs with equal size automatically satisfy momentum conservation on 
the macroscopic scale. Yet if we include explicitly the aqueous environment in which the ac- 
tomyosin network is immersed, hydrodynamic interactions between the nodes via the solvent 
should be taken into account. These interactions might modify the current simplified pic- 
ture and counteract any motor-induced force imbalance on individual nodes, thus validating 
momentum conservation on the microscopic scale. 

Dynamical evolution of the many-particle configuration {r^} due to these motor-driven 



nonequilibrium processes can be described by a master equation d^ /dt = I/^e^ ^^^ 
the configurational probability density ^({r^};t) with LNE^({^},i) = J ^idr[[K {{r^} -^ 
{r})^({r'}, t) — K{{r} -^ {^^D^d^}?^)] where the integral kernel K{{r\} -^ {^i}) encodes 
the probability of transitions between different crosslink/node configurations. Our earlier 
description^!^ of the motor kicking rate, fc, still applies to the current case of correlated 
kicks, i.e., 

k = K[e{AU) exp(-5^/3AC/) + e{-AU) exp{-SdPAU)], (1) 

where tv is the basal kicking rate and Su{sd) denotes motor susceptibility to energetically 
uphill (downhill) moves, except that the free energy change AC/ now arises from pairs of 
displacements. Explicitly we write 

xi^5{fi - Pi - kj)5{rj - Pj + kj)w [/(••• ,Pir" y^jr") -U{--- ,fi,-" ,fj,-") ^({Pj;t) 
-S{ri -r'i + lij)5{rj - r' j - lij)w [/(••• ,r;,--- , f^-, • • • ) - t^(- • • ^^'ir" ^^'j^'") *({^i};*)|- 

The factor 1/2 avoids double counting in the summation over all pairs. The quantity C^j, 
much like an element of a contact map in description of protein structures, defines whether 
the node pair {i^j) is connected by an active bond and thus subject to anti-correlated 
displacements (/^j, —lij)\ Cij = Cji = 1 for motor-bonded pairs while Cij = Cji = for non- 
bonded pairs. Our description of the rates gives w\Ui — Uf\ = Q{Uf — Ui) exp[— 5^t/3(C// — 
Ui)] + Q{Ui - Uf) exp[-Sd/3{Uf - Ui)]. 

Assuming symmetric motor susceptibility, i.e. Su = Sd = 5, one finds more simply 

i 3 
,,i^^Muin,r,)-Uir.-U,,rM.)\^^i^. . . , p. = f. - 4, ... , P • = f • + /"t., . . . }; t) 

_ -s^\u(ri+lij,rj-lij)-U(ri,f 



*({•••, rl,-.-,f„---};t)j. (2) 

We assume that kicks on different node pairs at any time are uncorrelated. The rates of 
possible kicking events depend on the instantaneous node configuration refiecting an assumed 
Markovian character of the dynamics. There is no angular average due to the definiteness 
of kicking directions for a given configuration. 



Note that the motor power strokes and thus kick steps are discrete occurring in a stochas- 
tic fashion. The correlated motions puU in slack locally while pulling taut neighboring fil- 
aments until a global balance is reached or a macroscopic collapse occurs, depending on 
whether the motors are downhill-prone (with a large positive s) or load-resisting (with a 
small or negative 5), respectively. The latter may be relevant to the contractile rachet-like 
behavior^^ that operates to incrementally drive cell shape change and deform tissues. 

B. Numerical translation: dynamic Monte Carlo simulation 

To realize the finite-jump Markov process described by the (chemical) master equation 
(Eq. [2j) we performed dynamic Monte Carlc^^ simulations on the model actomyosin network 
(Fig. [1]). In these simulations we generated initially a three-dimensional random lattice of 
volumeless nodes (mimicking the crosslinking proteins) and connected the nearest-neighbor 
nodes (defined by the first shell of the pair distribution function) with nonlinear elastic 
bonds22 (mimicking the act in filaments) at a given probability Pc. We then distributed the 
myosin motors uniformly to the bonds at a given probability Pa and obtained an active bond 
map. Considering the anti-correlated kicks along individual active bonds as chemical reac- 
tion channels, we adopted a stochastic simulation algorithm^^ to execute the moves following 
the stochastic process defined by the model motor kicking noise (Eq. [1]). For a sufficiently 
large system, kicking events on different node pairs are effectively decoupled (consistent with 
summing over independent reaction channels in the master equation) . Intermediate thermal 
moves between successive chemical moves obey Brownian dynamics^^ implemented via the 
position Langevin equation. 

The bond properties are given by the elasticity onset or relaxed length of actin filaments 
Lg = 1.2 and the effective stretch modulus /37 = 2. Since the relaxed length is larger than 
the mean node separation (set as the length unit), the initial homogeneous (but amorphous) 
network has a considerable fraction of ffoppy bonds. We assumed a relatively high basal 
kicking rate (/^ = 0.01) and a large kick step size (/ = 0.2) such that the dimensionless 
motor activity (defined as nt^ / D^ where D^ is the thermal diffusion constant) is close to 
1 and thus the strength of chemical noise is at least comparable to that of the thermal 
noise. The system size is A^ = 256 and periodic boundary conditions are applied. The 
relevant biophysical parameters to vary include network connectivity (Pc) as well as motor 



concentration (P^) and susceptibility {s). 

III. MAIN FINDINGS 

A. Role of motor susceptibility (s) and concentration (P^) in contractile behavior 

We first study how the concentration (P^) and susceptibility (s) of motors contribute to 
the collective behavior. As illustrated in Fig.O depending on the specific combinations of Pa 
and 5, distinct nonequilibrium structures emerge. For a force-percolating network (i.e. one 
with connectivity beyond the percolation threshold) when partially motorized (P^ < 1) un- 
der susceptible (large s) anti-correlated kicks, the active nodes (those with motor-attached 
bonds; shown as red spheres) begin to aggregate and tend to separate from the passive 
nodes (those with no active bonds; shown as blue spheres). The left panel of Fig. [2^ shows 
snapshots of both the initial and later node configurations. The corresponding developed 
network structure (Fig. [2^ right) exhibits clumps of fioppy bonds (concentrated short green 
lines) connected by tense bonds (long red lines). The overall rigidity (i.e. homogeneity 
on large scales) of the structure is protected by susceptible motors which tune the balance 
between local bond contraction and neighboring bond stretching such that energetically un- 
favorable tense states are avoided. At a low concentration of adamant motors (small 5), 
however, active nodes and their aggregates tend to "glue" together progressively the pas- 
sive nodes and their condensates (Fig. [2h). A finite spatial extent of the condensate and a 
non-vanishing fraction of taut bonds remain due to the insufficient cooperativity (low Pa) 
between local aggregation events. More dramatically, significant spatial heterogeneity forms 
when the system is driven by a large number of uphill-prone motors with negative suscepti- 
bility: the cooperative action of load-resisting motors induces a multistage aggregation and 
coarsening of the nodes (Fig. [2]3 upper row) finally leading to a macroscopic contraction 
of an initially homogeneous network into a dense clump of buckled filaments (Fig. [St lower 
row). This multistage coarsening process involves three steps: (1) local bond contraction and 
node aggregation giving fioppy clumps connected by tense filaments; (2) coarsening of the 
aggregates leading to filament alignment and formation of tense bundles; (3) coalescence of 
the larger aggregates into a single condensate accompanies collapse of the tense bundles into 
a fioppy clump. The aligned tense bundles formed before the eventual collapse constitute a 



taut state that can generate contractile forces. 

We will investigate the scenario of motor-driven aggregation and pursue its analogy to 
arrested phase separation later. In the present work we will focus on the regime for macro- 
scopic contraction. 

To identify the required motor properties for active contractility, we performed many 
simulations to obtain the evolution (in Monte Carlo time ^mc) of the statistical character- 
istics for a series of motor susceptibilities (Fig. Sla) and concentrations (Fig. Sib). These 
measures consistently indicate the existence of a threshold motor concentration ((Pa)th) 
and a threshold susceptibility (^th) for the onset of macroscopic contraction within the 
simulation time window. When s is less than 5th at an intermediate Pa value (Fig. Sla: 
s < 0.02, Pa = 0.5) or for Pa > {Pa)th at a smaU s value (Fig. Sib: Pa > 0.3,5 = 0), the 
fraction of taut bonds drops to essentially zero, indicating that an initially homogeneous 
percolating network collapses into a floppy clump. More interestingly, the total energy flrst 
rapidly rises and reaches a maximum before falling to zero with the fraction of taut bonds. 
This nonmonotonic behavior suggests that the system flrst works against an energy barrier 
due to the formation of a transient tense state (having highly stretched bundles induced 
by adamant or uphill-prone motor kicks) , and then cooperative action of suflSciently load- 
resisting motors drives the system over the barrier to allow energetically downhill moves 
via subsequent coarsening (as shown in Fig. [Jb). The rapid increase and saturation of the 
mean squared displacement (MSD) mirrors the evolution of the fraction of taut bonds and 
results from the formation of a single isolated floppy clump. Larger s and/or higher Pa 
yields a lower barrier to the collapsed state (as indicated by arrows in Fig. SI). Moreover, 
lower Pa necessitates having a smaller s to induce macroscopic contraction. In other words, 
weaker motor cooperativity requires a stronger load-resisting tendency to trigger contractile 
instability. When Pa < {Pa)th and/or s > 5th the structure remains homogeneous except for 
modest local node aggregation and network deformation. 

B. State diagram for active contractility: Interplay of network connectivity and 
motor cooperativity 

We now examine the interplay of network connectivity (Pc) and motor concentration (Pa) 
in forming contractile structures. The equivalence between the crosslink concentration and 
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the fraction of bonded neighboring crosshnks Pc (both are proportional to the number of 
inter-crosshnk segments and define the number of bond constraints), and between the motor 
concentration and the fraction of active bonds Pa (both are proportional to the number of 
crosslinks subject to motor kicks and determine the spatial cooper at ivity between motors) 
allows us to compare our state diagram constructed on the Pa-Pc plane with the experimental 
resullp^ shown for the parameter space of concentration ratios [myosin] /[actin] versus [a- 
actinin]/[actin]. 

We present our state diagram showing the dependence of macroscopic contractility on 
the network connectivity and motor concentration at a small s value {s = 0.01) in Fig. [31 
Red crosses denote contractile networks while blue circles denote non-contractile networks. 
By "contractile" we mean a complete collapse of an initially homogeneous network into a 
floppy clump within 10^ MC steps, monitored by the vanishing of total potential energy and 
the fraction of taut bonds. On top of the flgure, we display the initial network structures 
for several typical values of network connectivity Pc with the average number of bonded 
neighbors z. As observed experimentally,^^ we identify a threshold motor concentration and 
a window of network connectivity for active contractility, which deflne a parameter region 
as marked by the purple open frame. The new feature is a small window of motor concen- 
tration for contraction at high connectivity (marked by a green closed frame). These two 
aspects vividly demonstrate the interplay of network connectivity and motor concentration 
for global contraction: at any connectivity beyond the percolation threshold {Pc > 0.2), a 
sufliciently high motor concentration is required to achieve cooperativity among local con- 
traction events; on the other hand, since the bond constraints are strong at high connectivity, 
motor concentration cannot be too high since force asymmetry (imbalanced tug-of-war) is 
necessary to trigger local contraction. In the in vitro experiments,^^ no macroscopic contrac- 
tion was observed on the hour time scale at high crosslink concentrations. This observation 
is not incompatible with our results since inclusion of excluded volume eflFects should give a 
dramatic slow-down of the contractile dynamics owing to jamming and/or glass transition 
which would account for the absence of observable contraction on laboratory time scales. 

Outside the framed regions, at low motor concentration and/or low connectivity as well as 
at high motor concentration and high connectivity, there is no macroscopic contraction. We 
illustrate the failure of contractility at low connectivity (Fig. S2a) or low motor concentration 
(Fig. S2b). For Pc as low as 0.1, the average number of bonded neighbors is no greater than 
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1, lack of tension percolation thus prevents global contraction. This becomes more obvious 
as we increase the motor concentration; an increasing trend of local collapse with rising Pa 
value is apparent (Fig. S2a), leading to more compact aggregates and disconnected floppy 
clumps (several typical spots have been circled for the highest Pa case). At very low motor 
concentration P^ = 0.1 (Fig. S2b), formation of sparse and small active-node foci did not 
dramatically reshape the network, since rare and separate contraction events are insuflicient 
to trigger global contractile instability. 

A quantitative demonstration of the interplay between Pc and Pa is given in Fig. S3. 
Macroscopic contraction occurs either for high motor concentration (upper row, Pa = 0.7) 
at intermediate connectivity {Pc = 0.3,0.5) or for high connectivity (lower row, Pc = 0.7) 
at intermediate motor concentrations {Pa = 0.3,0.5). Lack of percolation at Pc = 0.1 is 
signalled by the diflFusive behavior of the mean square node displacement MSD linearly in- 
creasing with tMC (red line indicated by arrow in Fig. S3a right panel) since the disconnected 
aggregates merely undergo thermal motion. On the other hand, a small steady value for 
the MSD at P^ = 0.1 (red line indicated by arrow in Fig. S3b right panel) reflects the low 
cooperativity which causes no more than modest local distortions. A balanced tug-of-war 
at high Pc and high Pa disfavors even local deformations yielding a lower fraction of taut 
bonds, lower energy and smaller MSD as Pa increases. This can be seen by comparing the 
measures for P^ = 0.1 and Pa = 0.7 at Pc = 0.7 in Fig. S3b. 

C. Contrast with the case for uncorrelated kicks 

Our earlier study of a model cytoskeleton with uncorrelated isotropic kicks acting on in- 
dividual nodes^^ revealed that force percolation and mechano chemical coupling due to sus- 
ceptible motors can conspire to maintain a spontaneous flow, whereas adamant motor kicks 
promote fluidization (characterized by a vanishing localization strength of the nodes and 
formation of a disordered tense network). Under anti-correlated kicks acting along the lines 
of centers of motor-bonded node pairs, however, no vectorial flow transition is found even for 
a very large kick step size (/ c^ 0.5). This is probably because the restrictiveness of (local) 
kicking directions for a given conflguration impedes a global concerted net movement of the 
whole lattice. Instead, balanced local contraction and neighboring stretch in the presence 
of force percolation results in a network of floppy spots connected by tense fllaments which 
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remains homogeneous on large scales and bears some resemblance to arrested phase separa- 
tion occurring in low-packing-fraction physical gels with strong short-range attractions^^. At 
intermediate connectivity, force imbalance sensed by susceptible motors still induces phase 
separation into large floppy clumps connected by taut inter-clump bonds yet oscillations are 
no longer found. Fluidization is replaced by global contraction at above-threshold motor 
concentration. 

When the uncorrelated-kicking system was driven by uphill-prone motors with negative 
susceptibility, aster-like patterns formed (Fig. H] left). Instead, now, multistage coarsening 
and eventual macroscopic collapse occur (Fig. H] right), after surmounting a high energy bar- 
rier caused by transient tense states with stretched bundles. Comparing the bond tension 
patterns for these two cases in Fig.Hl we see that global contraction requires correlated move- 
ments that locally buckle a filament yet impose strain on the bonds at both ends. Without 
this correlation the aster pattern cannot collapse. We display in Fig. S4 the snapshots in the 
course of aster formation which clearly demonstrate that node aggregation and coarsening 
leads to progressively more tight filament bundling and bundle alignment, a mechanism that 
operates for both the correlated- and uncorrelated-kick cases. 

IV. CONCLUSION AND DISCUSSION 

We have simulated a microscopic model for the actomyosin cytoskeleton as a motorized 
cat's cradle that combines the asymmetric load response of individual actin filaments with 
correlated motor-driven events. This model reproduces the dependence of active contractility 
on microscopic parameters observed in reconstituted actomyosin networks. The simulations 
allow us to identify several necessary conditions for active contractility: highly asymmetric 
load responses of the filaments are needed for local contractile behavior under anti-correlated 
kicks, and a minimal network structure beyond percolation threshold is required to prop- 
agate local contraction. Sufliciently cooperative load-resisting motors manage to drive the 
system through energy-costly intermediate states and incorporate local buckling events into 
macroscopic contraction via a multistage coarsening process. 

This numerical study provides an explanation for the formation and contractile dynamics 
of disordered condensed state of actomyosin in vivo. The fact that such a simplifled model is 
able to mimic cellular self-organization states and their contractile dynamics suggests that 
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purely physical interactions contribute to the regulation of cell and tissue morphogenesis. 
Nevertheless specific biochemical signaling events certainly contribute to the localized as- 
sembly and activation of myosin foci such as occurs in cleavage furrow during cytokinesis 
and in wound borders. 

It is clearly necessary in the future to take into account the excluded volume effect which 
is expected to dramatically slow the contractile dynamics at high crosslink concentrations, 
thus accounting for the failure to see contractility on laboratory time scales. Long-range 
hydrodynamic interactions between the node pairs may also change the qualitative physics. 
Hydrodynamic correlations should facilitate propagation of nearby kicking events and break 
the continuous symmetry so as to allow spontaneous directed motion with motor-driven 
hydrodynamics-mediated pulsed contractions. Finally, in vivo one must incorporate motor 
attachment and detachment and crosslink binding and unbinding for complete realism. 

Support from the Center for Theoretical Biological Physics sponsored by the National 
Science Foundation (Grant PHY-0822283) is gratefully acknowledged. 
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FIG. 1: Schematic of the model system: A cat's cradle composed of a three dimensional amor- 
phous network of crosslinked nonlinear-elastic bonds (purple lines) where motor-driven anti- 
correlated kicks induce pulsed local bond contraction and node aggregation. Spheres represent 
nodes/crosslinks and the yellow fuzzy objects stand for myosin II motor proteins. The size of 
the nodes and motors is exaggerated; excluded volume interaction is not implemented in current 
simulations. Red nodes are active nodes having motor-attached bonds. These are subject to motor 
kicks. Blue nodes are passive nodes only undergoing Brownian motion. An enlarged view of a unit 
of local contraction (circled) shows a pair of motor-bonded crosslinks undergoing anti-correlated 
kicks with fixed step size / along the line connecting their centers where fij is a unit vector pointing 
from node i to node j. 
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Multistage coarsening/aggregation 
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FIG. 2: Illustration of how motor concentration Pa and susceptibility 5 contribute to their col- 
lective behavior, (a) Left: motor-driven aggregation of active nodes (red spheres) among the 
passive nodes (blue spheres); right: formation of dense floppy clumps (concentrated short green 
lines) inter-connected by tense bonds (long red lines). Circles mark the corresponding regions 
of node aggregation and bond collapse. Pc = 0.2, P^ = 0.5,5 = 1. (b) Active nodes tend to 
"glue" together passive nodes and their aggregates. Pc — 0.5, P^ — 0.2,5 = 0. (c) Multistage 
coarsening/aggregation of active condensates driven by high-concentration (P^ = 1) uphill-prone 
{s = —0.5) motors. Pc = 0.5. Shown are the temporal evolution of the node conflguration (upper 
row) and of the corresponding network structure (lower row). 
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FIG. 3: State diagram showing the dependence of active contractihty on motor concentration Pa 
and network connectivity Pc- 5 = 0.01. Red crosses denote contractile networks and blue circles 
denote non-contractile networks. Top row displays the initial homogeneous network structures with 
increasing connectivity Pc and average number of bonded neighbors z. Macroscopic contraction 
occurs in the two framed regions: (1) purple open frame: intermediate Pc and above-threshold Pa] 
(2) green closed frame: high Pc and intermediate Pa- 
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Uncorrelated kicks 



Correlated kicks 



FIG. 4: Aster versus condensate. Pc = 0.5, P^ = 1,5 = —0.5. Left: Under uncorrelated kicks, 
an initially disordered and homogeneous network self-organizes into highly tense and ordered 
"asters", composed of tense bundles radiating from the junctions where floppy bonds concentrate 
(see the zoom-in image). Right: Driven by anti-correlated kicks, an initially homogeneous force- 
percolating network flrst develops transient tense states consisting of highly-stretched bundles, but 
then abruptly collapses into a single floppy clump. 
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Threshold motor susceptibility (s) for contractility 
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Fig. SI I Thresholds for active contractihty. (a) Evolution of the statistical measures 
for a series of motor susceptibilities {s = 0,0.01,0.02,0.05 and 0.1). Pc = 0.5, P^ = 0.5. 
When s < 0.02 the fraction of taut bonds and the total energy drop to essentially zero 
after surmounting an energy barrier due to tense intermediates. Mean square displacement 
(MSD) mounts to a plateau as a consequence of the formation of a single floppy clump. 
Smaller s yields a higher barrier and faster collapse (indicated by arrows), (b) Evolution of 
the measures for a series of motor concentrations {Pa = 0.1,0.2,0.3,0.4,0.5 and 1). Pc = 
0.5, s = 0. When Pa > 0.3 global contraction occurs. Larger Pa leads to lower barrier and 
faster collapse (indicated by arrows). 
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Fig. S2 I Illustrations of the failure of contractility at low network connectivity 
or low motor concentration. 5 = 0. (a) Pc = 0.1. Lack of tension percolation due 
to low degree of bonding prevents global contraction. As Pa increases (left to right: Pa = 
0.1, 0.5, 0.8 and 1) system exhibits increasing trend of local aggregation (lower row) resulting 
in an increased fraction of floppy bonds (upper row). When driven by high-concentration 
(Pa = 1) adamant motors, the initially homogeneous network exhibits local collapses into 
disconnected clusters of buckled filaments (circles mark the typical regions), (b) Pa = 
0.1, Pc = 0.5. Sparse and modest local network distortion (upper) and node aggregation 
(lower) are insuflicient to trigger global contractile instability. 
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Interplay of Pc and Pa in macroscopic contraction 
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Fig. S3 I Interplay of network connectivity and motor concentration in macro- 
scopic contraction, s = 0.01. (a) Pa = 0.7, Pc = 0.1,0.3,0.5,0.7: contraction occurs for 
intermediate Pc. Larger Pc yields higher "barrier" to coUapsed state since more bond con- 
straints results in more intense tug-of-war. (b) Pc = 0.7, Pa = 0.1,0.3,0.5,0.7: contraction 
occurs for intermediate Pa- Larger Pa lowers the barrier since cooperativity between local 
contraction events is enhanced. 
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Aster formation 




Fig. S4 I Snapshots in the course of aster formation. Upper row shows the bond 
structure: filament bundhng evolves from loose to tight. Lower row presents the corre- 
sponding node configuration: aggregation toward the corners of the simulation box becomes 
progressively more compact. Pc = 0.3, P^ = 1,5 = —0.5. 
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